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The results of a joint NASA/Army/Bell Helicopter Textron wind-tunnel test to assess the potential of Gener- 
alized Predictive Control (GPC) for actively controlling the swashplate of tiltrotor aircraft to enhance aeroelastic 
stability in the airplane mode of flight are presented. GPC is an adaptive time-domain predictive control method 
that uses a linear difference equation to describe the input-output relationship of the system and to design the con- 
troller. The test was conducted in the Langley Transonic Dynamics Tunnel using an unpowered 1 /5-scale semispan 
aeroelastic model of the V-22 that was modified to incorporate a GPC-based multi-input multi-output control algo- 
rithm to individually control each of the three swashplate actuators. Wing responses were used for feedback. The 
GPC-based control system was highly effective in increasing the stability of the critical wing mode for all ot the 
conditions tested, without measurable degradation of the damping in the other modes. The algorithm was also ro- 
bust with respect to its performance in adjusting to rapid changes in both the rotor speed and the tunnel airspeed. 
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Notation 

matrices (formed from a, and p,) in multi- 

step output prediction equation 

control horizon 

prediction horizon 

objective function to be minimized 

time index 

number of data time steps used for SID 
number of feedback outputs to GPC 
order of ARX equation describing system 
weighting matrices for inputs and outputs 
number of control inputs from GPC 
vector of digitized input time histories 
vector of computed control inputs 
matrix of input and output measurements 
weighting factors for inputs and outputs 
matrix of OMP determined by SID 
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y(k) vector of digitized response time histories 

a l \ P control law gain matrices 

a„ P coefficient matrices (OMP) appearing in 

ARX equation and determined by SID 
Si blade kinematic pitch-flap coupling angle 

e error function (difference between target 

and predicted responses) 

Abbreviations: 

ARX autoregressive with exogenous input 

MIMO multi-input/multi-output 

OMP observer Markov parameters 

SID system identification 

INTRODUCTION 

Tiltrotor aircraft operating at high speeds in the 
airplane mode of flight are susceptible to a propro- 
tor/pylon instability akin to propeller whirl flutter. Such 
an instability was first encountered during full-scale 
testing of the Bell XV-3 tiltrotor in the NASA-Ames 
40- by 80-foot Wind Tunnel in 1962. Following this 
incident, extensive analytical and experimental studies 
of small dynamically-scaled models were initiated by 
Bell with the threefold objective of providing a physical 
understanding of the phenomenon, developing an ana- 
lytical method to predict such instabilities, and identify- 
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ing corrective design changes. By 1965, both an expla- 
nation for and a means for eliminating the instability 
were found (refs. 1-2). The XV-3 aircraft was then 
modified to incorporate changes indicated by analyses 
and retested successfully in the NASA-Ames 40- by 80- 
foot Wind Tunnel in 1966. The proprotor/pylon dy- 
namic studies initiated by Bell in support of the XV-3 
investigation were followed by several other studies 
conducted by government and industry over the next 
decade (see, for example, refs. 3-6). Taken as a whole, 
these studies helped to establish the technology base 
needed to later successfully address the issue of propro- 
tor/pylon/wing aeroelastic stability in the design of the 
XV- 15 tiltrotor research aircraft in the mid-1970s and 
the V-22 Osprey in the mid- 1980s. 

The mechanism of the proprotor/pylon instability 
experienced by the XV-3 was identified by Bell during 
the investigation of that incident and was reported by 
Hail (ref. 1). The primary destabilizing factor was 
found to be hub inplane shear forces resulting from ro- 
tor precession when operating at high forward speeds in 
the airplane mode of flight. The hub shears identified as 
being responsible for the subject instability arc actually 
a subset of a larger set of hub forces and moments that 
act on a tiltrotor aircraft (ref. 4). Disturbances occurring 
in flight can excite either the elastic or rigid-body 
modes of an aircraft in an oscillatory manner. For a 
tiltrotor aircraft, any motions of this type effectively 
represent oscillatory translational and rotational motions 
of the proprotor shaft in space. This leads to proprotor- 
generated aerodynamic forces and moments that are a 
function of these oscillatory motions. Figure 1, from 
reference 4, shows the perturbation rotor-induced aero- 
dynamic forces acting at the hubs of a tiltrotor aircraft 
executing small pitching and yawing motions when op- 
erating in an airplane mode of flight. From the position 
of these forces on the aircraft, it is clear that the forces 
shown can influence aircraft longitudinal and lateral- 
directional stability. However, the shear forces H and Y 
can, quite independently of any aircraft rigid-body mo- 
tions, also destabilize the proprotor-pylon-wing system 
aeroelastically. These are the shear forces identified by 
Hall (ref. 1) as the drivers for proprotor-pylon instabil- 
ity. They are a direct consequence of airload moments 
that are generated to precess the rotor in space in re- 
sponse to shaft pitching and yawing motions. In addi- 
tion to a true whirl instability involving both pitching 
and yawing motions of the pylon, proprotor/pylon insta- 
bility can occur in a single plane (either pitch or yaw), 
depending on the pylon support stiffnesses. Propel- 
ler/pylon whirl flutter, on the other hand, is driven by 
aerodynamic cross-stiffness moments and can occur 
only if there is pylon flexibility in both pitch and yaw. 
A discussion of these and other important differences in 
the aeromechanical behavior of propellers and propro- 
tors is given in reference 4. 



Figure 1. - Perturbation rotor-induced aerodynamic 
forces acting on a tiltrotor aircraft during pitching and 
yawing oscillations (rotors interconnected). 

Proprotor/pylon aeroelastic stability continues to 
be a primary consideration in the design of tiltrotor air- 
craft. Because wing torsional stiffness is the major 
structural design parameter influencing this phenome- 
non, the stability requirements of current tiltrotor air- 
craft (XV- 15, V-22, BA-609) have been attained by 
using thick, torsionally stiff wings having a 23% thick- 
ness-to-chord ratio. Such wings provide the torsional 
stiffness required for stability in those aircraft, albeit at 
the expense of cruise efficiency and maximum speed. 
The use of thinner wings would permit higher cruise 
speeds, increased range, and improved productivity. 
However, the attendant reduction in wing stiffness 
would bring with it the problem of proprotor/pylon in- 
stability. Thus, aeroelastic instability of the propro- 
tor/pylon/wing system stands as a major barrier to in- 
creasing the maximum speed capability of tiltrotor air- 
craft. Both passive and active methods for extending 
proprotor/pylon stability boundaries for conventional 
tiltrotor aircraft have been investigated. References 7- 
12, which describe studies into the use of wing and 
blade aeroelastic tailoring techniques for improving 
stability, are illustrative of passive approaches. 

While there have been a number of studies deal- 
ing with the use of active controls for improving the 
aeroelastic behavior of tiltrotor aircraft, most of these 
have addressed the problem of gust and maneuver load 
alleviation and there has been only limited attention 
given to the use of active controls for stability augmen- 
tation (refs. 13-16). Reference 13 investigated the ap- 
plication of swashplate feedback for augmentation of 
aeroelastic stability as part of a broader study of feed- 
back control for improving the aeroelastic and rigid- 
body flight characteristics of tiltrotor aircraft. The base- 
line configuration of those studies was a Boeing soft- 
inplane tiltrotor known as the Model 222. Bode analy- 
ses were used to define the appropriate gains and phases 
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to be applied to the wing responses that were to be fed 
back to the swashplate cyclic inputs. References 14-15 
were analytical studies of feedback control for increas- 
ing proprotor/pylon stability on an XV- 15 size tiltrotor 
aircraft. A feedback gain matrix was introduced into the 
formulation by adding a feedback loop to the equations 
of motion linearized about a flight condition of interest. 
Wing tip vertical velocities and accelerations were used 
for feedback. The gains needed to stabilize the system 
were determined by simply varying the terms in the gain 
matrix until an eigenvalue analysis of the closed-loop 
system indicated a stable system. Reference 16 was an 
analytical study into the use of linear quadratic regulator 
(LQR) techniques for determining the wing feedback 
gains needed to stabilize the whirl modes of a tiltrotor 
aircraft using an active swashplate. The method was 
studied using a mathematical model that had been de- 
veloped earlier for a full-size semi-span configuration of 
the XV- 15. Control design was done in modal space. 
The discrete state-space equations were transformed to 
modal form to allow for separation of stable and unsta- 
ble modes during the design of the controller. A Kal- 
man-Bucy filter was employed as the state estimator to 
account for disturbances and noise. 

The Aeroelasticity Branch (AB) and the Army 
Research Laboratory's Vehicle Technology Directorate 
at NASA Langley Research Center, in collaboration 
with Bell Helicopter Textron Inc (BHTI), is evaluating 
an adaptive control technique known as Generalized 
Predictive Control (GPC) (ref. 17) to assess its potential 
for actively controlling the swashplate of tiltrotor air- 
craft to enhance aeroelastic stability in both helicopter 
and airplane modes of flight. The term adaptive as used 
here refers to a control system in which the parameters 
describing the model are identified under operational 
conditions using measured input/output data, and the 
updated system parameters are then used to compute a 
new set of control gain matrices and the next set of 
commands to be sent to the swashplate actuators, with 
all computations being done on-line. GPC is a time- 
domain method that uses a linear difference equation to 
describe the input-output relationship of the system and 
to design the controller. An ARX-type difference equa- 
tion is being employed in the present study. The ARX 
model is used for both system identification and control 
design. The coefficient matrices of the ARX equation 
are the quantities determined by the identification algo- 
rithm. Closed-loop control is enhanced by performing 
the system identification in the presence of the external 
disturbances acting on the system. The coefficients of 
the ARX model are assembled into a multi-step (finite- 
horizon) output prediction equation, the desired (target) 
response is specified, and the resulting expression is 
used to form an objective function. Minimization of the 
objective function leads to an expression for the control 
to be applied to the system. 


A 1 /5-scale, semispan aeroelastic model of the 
V-22 designed and built by Bell in 1981 is being used in 
the experimental evaluation of GPC. That model, on 
loan to the AB from the Navy, has been refurbished to 
form a tiltrotor research testbed called the Wing and 
Rotor Aeroelastic Test System (WRATS) lor use in the 
Langley Transonic Dynamics Tunnel (TDT) (fig. 2). 



Figure 2. - WRATS tiltrotor testbed installed in TDT. 


Three exploratory experimental studies to investigate 
the use of a multi-input/multi-output (MIMO) GPC- 
based method for active control have been conducted on 
the WRATS model to-date (ref. 18). These include a 
hover test with the baseline stiff-inplane gimbaled rotor 
(May 1998), a ground resonance test employing a soft- 
inplane variant of the gimbaled rotor (October 1999), 
and a wind-tunnel test with the baseline stiffness rotor 
(April 2000). 

The purpose of this paper is to summarize the re- 
sults of the joint NASA/Army/BHTI investigation that 
was conducted in the TDT in April 2000 in the initial 
assessment of the potential of GPC for augmenting sta- 
bility in tiltrotor aircraft operating in the airplane mode 
of flight. The WRATS model was modified to incorpo- 
rate a GPC system for independently controlling the 
actuators of the fixed-system swashplate. Active con- 
trol was introduced into the fixed-system swashplate 
using three high-frequency servo-controlled hydraulic 
actuators mounted aft of the swashplate inside the pylon 
fairing. The actuator commands included the steady 
pilot commands as well as the actuator motions called 
for by the active control algorithm. Wing bending and 
torsion strain gages located near the wing root were the 
sensors used to provide the response measurements 
needed by the active control algorithm. The GPC- 
computed commands were summed with the pilot trim 
commands to produce the desired swashplate actuation 
about the trim state. The effectiveness of the rotor 
swashplate in increasing stability (damping) was inves- 
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tigated over a range of tunnel airspeeds for a single rotor 
rotational speed. The robustness of the active control 
system to rapid variations in tunnel airspeed and rotor 
speed, both singly and in combination, was studied. 
The effect of GPC on blade loads and pitch link loads 
was not assessed. 

The pertinent equations underlying the method 
are presented and discussed and the strategy of com- 
puter implementation is described, including system 
identification, calculation of control law matrices, and 
calculation of the control commands sent to the swash- 
plate actuators. Considerations related to implementa- 
tion for on-line computations are discussed. Experi- 
mental results are then presented illustrating the stability 
augmentation that was obtained with the GPC-based 
active control system. 


GENERALIZED PREDICTIVE CONTROL 
Introductory Remarks 

Predictive controllers were initially introduced in 
the chemical industries for controlling chemical proc- 
esses and have found applications in a wide variety of 
industrial processes (e.g., ref. 19). Predictive control 
refers to a strategy wherein the decision for the current 
control action is based on minimization of a quadratic 
objective function that involves a prediction of the sys- 
tem response at some number of time steps into the fu- 
ture. A variety of predictive controllers have been pro- 
posed (e.g., ref. 20). Among these. Generalized Predic- 
tive Control (GPC), which was introduced in 1987 (ref. 
17), has received notable attention by researchers. GPC 
is a time-domain multi-input-multi-output (MIMO) pre- 
dictive control method that uses a linear difference 
equation to describe the input-output relationship of the 
system. The input-output equation is used to form a 
multi-step output prediction equation over a finite pre- 
diction horizon while subject to controls imposed over a 
finite control horizon. The control to be imposed at the 
next time step is determined by minimizing the devia- 
tion of the predicted controlled plant outputs from the 
desired (or target) outputs, subject to a penalty on con- 
trol effort. 

A novel version of the GPC procedure was de- 
veloped at NASA Langley Research Center in 1997 for 
efficient computation and unknown disturbance rejec- 
tion by Dr. Jer-Nan Juang and his co-workers. Their 
work has resulted in a suite of MATLAB m-liles that 
have been collected into a Predictive Toolbox that can 
be used by researchers for GPC studies. A summary of 
the SID and control theory underlying their develop- 
ment is found in references 21-29, among others. 


The essential features of the adaptive control 
process used in the present GPC investigation are de- 
picted in figure 3. The system has r control inputs w, m 
measured outputs y, and is subject to unknown external 



Figure 3. - Block diagram of identification and control 
procedure. 


disturbances d. Measurement noise is also present. 
There are two fundamental steps involved: (1) identifi- 
cation of the system, and (2) use of the identified model 
to design a controller. As mentioned earlier, a linear 
difference equation or model is used to describe the 
relationship between the input and the output of the sys- 
tem. A linear input-output model gives the current out- 
puts as a linear combination of past input and output 
measurements. The input-output model used in the pre- 
sent work has the form of what is referred to as an ARX 
(autoregressive with exogenous input) model. The 
ARX model is used for both system identification and 
controller design. System identification is done on-line, 
but not at every time point, in the presence of any dis- 
turbances acting on the system. In this way, an estimate 
of the disturbance model is reflected in the identified 
system model and the disturbance does not have to be 
modeled separately. This approach represents a case of 
feedback with embedded feedforward. Because the 
disturbance information is embedded in the feedforward 
control parameters, there is no need for measurement of 
the disturbance signal (ref. 28). The parameters of the 
identified model are used to compute the predictive con- 
trol law. A random excitation u i( j (sometimes called 
dither) is applied initially with u ( . equal to zero to iden- 
tify the open-loop system. Dither is added to the 
closed-loop control input u c if it is necessary to re- 
identify the system while operating in the closed-loop 
mode, 

Form of Model and Control Law Equations 

The relationship between the input and output 
digitized time histories of a MIMO system are described 
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where 


by an ARX model that has the form 


y(k) = a.y(k-\) + a^y(k- 2) + ■•• +a y{k-p) 

(I) 

+ fi f) u(k) +fi ] u{k-\) + ■■■ +p l ,u(k-p) 

This equation states that the current output y(k) at time 
step k may be estimated by using p sets of the previous 
output and input measurements, y(k- 1 y(k-p) and 
u(k-I ),...,u(k-p)\ and the current input measurement 
u(k). The integer p is called the order of the ARX 
model. The coefficient matrices a { and appearing in 
this equation are referred to as observer Markov pa- 
rameters (OMP) or ARX parameters and are the quanti- 
ties to be determined by the identification algorithm. 
Closed-loop robustness is enhanced by performing the 
system identification in the presence of the external 
disturbances acting on the system, thereby ensuring that 
disturbance information will be incorporated into the 
system model. The goal of SID is to determine the 
OMP based on input and output data. The OMP may be 
determined by any SID technique that returns an ARX 
model of the system. 

The ARX model is used to design the controller 
and leads to a control law that in the case of a regulator 
problem has the general form given by 

U (k) = a{ y(k-l) + a[ y(k-2) + ■ • * +a L „ y( k-p) 

1 (2) 

+/?; u(k- 1) +/?2 y(k- 2) + • • • + P],u{k-p) 

Equation 2 indicates that the current control input u c (k) 
may be computed using p sets of the previous input and 
output measurements. The coefficient matrices a and 
ft appearing in this equation are the control gain matri- 
ces. 


System Identification 

System identification in the presence of the op- 
erational disturbances acting on the system is the first of 
the two major computational steps. The external distur- 
bances acting on the system are assumed to be unknown 
(unmeasurable). The number of control inputs is r and 
the number of measured outputs is m. The system is 
excited with band-limited white noise for SID. These 
random excitations are input to all r control inputs si- 
multaneously and the m responses are measured. The 
digitized input and output time histories ( u and y) at / 
time points are then used to form the data matrices y and 
V in the equation 


y = r v/ 


(3) 


y- 

,v(0) 

V(l) V(2) •• 

■y(p) ■■ 

■yd- 1) 


(4) 
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v(p-\) 

- v(/ 

-2) 


V= 


v(0) 

v(p-2) 

v(/ 

-3) 

(5) 




v(0) 

••• v(/- 

P~ 1)| 



[r + (r + m)p] x / 


Equations 4 and 5 follow from writing the discrete-time 
state-space equations for a linear time-invariant system 
at a sequence of time steps k = 0, 1,..., M and grouping 
them into matrix form. The vector v(k) appearing in the 
data matrix V is formed from the vectors u(k) and y(k) 
according to 


v(0= 


u(k) 

y(k) 


( r + m ) x 1 


( 6 ) 


The order of the ARX model (p) and the number of time 
steps (/) must be specified by the user. Some guidelines 
for their selection are given later. The sizes of the vec- 
tors and arrays are noted. 


In forming the matrices given in equations 4 and 
5, it has been assumed that the state matrix A is asymp- 
totically stable so that for some sufficiently large p , A k = 
() for all time steps k > p, and that an observer has been 
added to the system. It is through these expedients that 
the matrix V is reduced to a size amenable for practical 
numerical computation of its pseudo-inverse. The SID 
process yields OMP rather than system Markov parame- 
ters (SMP) because of the inclusion of an observer. A 
complete discussion of these aspects of the development 
may be found in reference 2 1 . 


Y is the matrix of observer Markov parameters 
that are to be identified and has the form 


y= 


/?0 $2 ^3 ftp 


mx r mxr m x m m xr mx tn m xr mx m m x r m x m 
;»x[/+ (r + /fj)/>} 


(7) 
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The solution for Y is obtained by solving equation 3 
for Y according to 

Y = y V f =y V T [V V T J' (8) 

where j denotes the pseudo-inverse. If the product V V T 
is a well-conditioned matrix of reasonable size, the or- 
dinary inverse can be taken as shown. Otherwise, a 
pseudo-inverse must be used. It should be noted that 
because the size of V V T is much smaller than V , a 
pseudo-inverse might be appropriate even if the product 
is well conditioned. 

Multi-Step Output Prediction Equation 

The one-step ahead output prediction equation given 
in equation 1 is the starting point for deriving the multi- 
step output prediction equation that is needed for de- 
signing a GPC controller. Using equation 1, the output 
at time step k+j may be written in the form 

y(k + j) = a\ J) y(k-l) + a\ j) y{k-2) + •• +a { * ) y(k-p) 

+ p {) u(k + j)+P!''u(k + j- 1) +-- + #VA) (9) 
+ Pi 1 hi (k - 1 )+ pfu (k- 2 ) 4 - ••• +/?*/ hi (k -p ) 

where the coefficient matrices are given by recursive 
expressions involving the matrices a t and appearing 
in the ARX equation (ref. 25). The system identifica- 
tion process described earlier determines the matrices 


and Pi. Equation 9 shows that the output y(k+j) at time 
step k+j may be estimated by using p sets of the previ- 
ous output and input measurements, y(k-\ y(k-p) 

and u(k-\) u(k-p), and the (unknown) current and 

future inputs, u(k), u(k+\),...,u(k+j). The GPC algo- 
rithm is based on system output predictions over a finite 
horizon h p known as the prediction horizon. To predict 
future plant outputs, some assumption needs to be made 
about future control inputs. In determining the future 
control inputs for GPC, it is assumed that control is ap- 
plied over a finite horizon h c known as the control hori- 
zon. Beyond the control horizon the control input is 
assumed to be zero. In GPC, the control horizon is al- 
ways equal to or less than the prediction horizon. Let- 
ting) in equation 9 range over the set of values j - 1,2, 
/? /r l, the resulting equations can be assembled into a 
multi-step output prediction equation having the form 



GO) 


The coefficient matrices T, B , and A are formed from 
combinations of the observer Markov parameters a, and 
Pi. The quantity \i, r (k) is the vector containing the fu- 
ture outputs, whereas u fw (k) is the vector containing the 
future control inputs yet to be determined. The quanti- 
ties u r (k-p) and y r (k-p) are vectors containing the previ- 
ous p sets of control inputs and outputs, respectively. 
The expanded form of this multi-step output prediction 
equation is shown in equation 1 1 . 
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The OMP a, and /?, determined in SID form the first 
block rows of the coefficient matrices Z A , and B in 
equation 1 1. The terms in the remaining rows are com- 
puted using the recursive relations indicated in the 
boxes (ref. 25). All terms in equation 1 1 are known, 
except for the h L sets of future commands and the h p sets 
of predicted responses. The goal of the GPC control 
algorithm is to determine the set of future commands 

u(k\ u(k+\) u{k+h c -\) that are required to achieve 

a desired predicted response y(k), y(k+ 1 ), .... y(k+h v - 1 ). 

It should be remarked that the system Markov 
parameters (SMP), which are commonly used as the 
basis for identifying discrete-time models for linear dy- 
namical systems, form the first block column in the ma- 
trix Z the remaining block columns are formed from 

subsets of the SMP. The Markov parameters are the 
pulse response of a system and are unique for a given 
system. The discrete-time state-space matrices A . B, C, 
and D are embedded in the SMP. 

Derivation of Control Law 

The predictive control law is obtained by mini- 
mizing the deviation of the predicted controlled re- 
sponse (as computed from the multi-step output predic- 
tion equation) from a specified target response over a 
prediction horizon h f) . To this end, one first defines an 
error function that is the difference between the desired 
(target) response y r (k) and the predicted response 
)'h P (k): 


u h (k) -~{T r RT +Q)'T I R [-y T (k)+Bu,(k-p)+Ay p (k-p)) 
= 7 ( -y T (k )+Bu r (k- p)+Ay p (k-p)) 


as the control sequence to be applied to the system over 
the next h c time steps. However, only the first r values 
(corresponding to the first future time step) 


u c ( k ) = -f y T ( k ) + P c u p (k -p) + a c y p ( k - p ) 

r x 1 

(15) 

arc applied to the r control inputs, the remainder are 
discarded, and a new control sequence is calculated at 
the next time step. 

The target response is zero for a regulator prob- 
lem and non-zero for a tracking problem. The matrix Q 
must be tuned to ensure a stable closed-loop system. 
Typically, h c is chosen equal to h p . However, h c may be 
chosen less than h r resulting in a more stable, but slug- 
gish, regulator. 

An expression for estimating the order of the 
ARX model that is to be used for SID is given by 


p>ceil 


number of 
sy stem states 


number of 
disturbance states 
m 


(16) 


j 


£ = y T (k)-y llp (k) 

= y r (k)-Tu, k ( k)-J3u p (k - p)-Ay p (k-p) 

( 12 ) 

An objective function J quadratic in the error and the 
unknown future controls is then formed: 

J = £ T R£ + ul Q U, k (13) 

Two weighting matrices are included in the objective 
function: Q (symmetric and positive definite) is used to 
weight the control effort and stabilize the closed-loop 
system; R (symmetric and positive semi-definite) is 
used to weight the relative importance of the differences 
between the target and predicted responses. Typically, 
Q and R are assumed to be diagonal and for Q to have 
the same value vv t . along its diagonal and R to have the 
same value vvy along its diagonal. Minimizing J with 
respect to u h Jk) and then solving for u hc (k) gives 


where ceil denotes rounding up the value of the quantity 
in parentheses to the next higher integer. The number 
of system states is typically chosen to be twice the num- 
ber of significant structural modes; the number of dis- 
turbance states is set to twice the number of frequencies 
in the disturbance; m is the number of output measure- 
ments. If measurement noise is of concern, the order of 
p so computed should be increased to allow for compu- 
tational poles and zeros to improve system identification 
in the presence of noise. In practice, simply choosing p 
to be 5-6 times the number of significant modes in the 
system is often adequate. The prediction and control 
horizons are set according to the relations 

h p > p h c <h p (17) 

Although h p can be set equal to p , h p is typically set 
greater than p to weight the control effort so as not to 
saturate the control actuators. If the control horizon is 
greater than the system order a minimum energy (mini- 
mum norm) solution is obtained wherein the com- 
manded output is shared so that the control actuators 
don’t fight each other. If h p is set equal to p one obtains 
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a so-called deadbeat controller (ref. 25). By extending 
the prediction and control horizons (and hence p , the 
order of the ARX model) to very large values, the GPC 
solution approaches that of the linear quadratic regulator 
(LQR). Thus, GPC approximates an optimal controller 
for large p (ref. 17). 

Weighting matrices Q and R are used to weight 
the control effort and to weight the relative importance 
of the differences between the target and predicted re- 
sponses, respectively. As mentioned earlier, Q and R 
are usually assumed to be diagonal and for Q to have 
the same value w c along its diagonal and R to have the 
same value w r along its diagonal. The control weight 
must be tuned to produce an acceptable solution without 
going unstable. Reducing w v increases control authority 
and performance but eventually drives the system un- 
stable. 

Key Features of Present Method 

An ARX model is employed to represent the sys- 
tem and is used for both system identification and con- 
troller design. The SID process used makes recourse to 
an observer to enable numerical computation of the 
pseudo-inverse needed for calculation of the OMP that 
comprise the coefficients of the ARX equation. The 
controller is thus inherently observer-based but no ex- 
plicit consideration of the observer needs to be taken 
into account in the implementation. In practice, the 
disturbances acting on the system are unknown or un- 
measurable. However, as discussed in reference 28, by 
performing the SID in the presence of the external dis- 
turbances acting on a system, a disturbance model is 
implicitly incorporated into the identified observer 
Markov parameters. However, the identified model 
must be larger than the true system model to accommo- 
date the unknown disturbances. Thus, the effects of the 
unknown disturbances acting on the system are embed- 
ded in the matrices A and B , and hence the control law 
matrices a and ff . If only the disturbances acting on 
the system change, there is no need to recalculate 
T (and hence / ) because T is formed solely from the 
SMP, which are unique for a given system (ref. 28). 
The solution for Y indicated in equation 8 involves 
forming the matrix products yV and VV . Here, these 
products are obtained using the computationally effi- 
cient procedure described in reference 23. 

Computational Considerations 

Several considerations dealing with computa- 
tions should be kept in mind when developing algo- 
rithms for GPC applications. 


SID should be done with the external distur- 
bances acting on the system so that information about 
the disturbances is embedded in the OMP. However, 
depending on the nature of the external disturbances, it 
may be possible to perform a SID on a system without 
the external disturbances and still determine a control 
law that results in satisfactory closed-loop performance. 

The computation of pseudo-inverses should be 
performed using Singular Value Decomposition (SVD) 
because of the latter’s ability to deal with matrices that 
are numerically ill conditioned. The use of pseudo- 
inverses (via SVD) is recommended even in cases 
where the ordinary inverse may seem appropriate (such 
as in the operation (VV 7 ) 1 indicated in equation 8). 

A number of decisions also need to be made. 
For example, should the system be re-identified and the 
control law matrices recalculated in every sampling 
period, every specified number of time steps (block up- 
dating), or only if some event requiring a re-ID occurs? 
Should all the calculations be done in real time, or can 
some be done in near real time? Should updating of the 
OMP be done in batch mode or recursively? 

Microprocessor speeds are such that it is now of- 
ten possible to complete the full cycle of GPC computa- 
tions and apply the commands to the actuators within 
one sampling period. However, the need for efficient 
computational algorithms and attendant coding is not 
expected to diminish. 

Implementation Considerations 

Several considerations must be taken into ac- 
count when actually implementing GPC algorithms in 
hardware for active controls work. 

The measured response time histories must be 
passed through a low-pass filter with a cut-off frequency 
set equal to the Nyquist frequency f N . The latter is cho- 
sen so that the maximum frequency of interest is about 
75% of/v The sampling frequency f s should be at least 
twice //y to prevent aliasing. However, if fs is made too 
large the low frequency modes will be poorly identified 
due to a loss of frequency resolution. A sampling rate 
between 2 to 3 times f N is generally sufficient. Once the 
sampling frequency has been selected, the minimum 
number of data points that should be used for SID fol- 
lows from the requirement of having 5-10 cycles of the 
lowest frequency mode in the measured response time 
histories. Normalization of the input and output data 
that is used for SID on the maximum actual or expected 
values of the data is often helpful numerically. The 
procedure will depend on whether the computations are 
being done in batch mode or recursively. 
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The computing tasks can be distributed among 
computers or different CPUs on a single computer. The 
choice will influence the extent of user involvement. 
The values of the input and output data that are being 
used during closed-loop operations must be carefully 
monitored to ensure that they fall within acceptable 
bounds. This is easily done using IF/THEN-type checks 
in the code. 

EXPERIMENTAL SETUP 
Wind Tunnel 

The experiment was conducted in the Langley 
Transonic Dynamics Tunnel, which is a continuous 
flow, single return, variable pressure tunnel having a 
test section 16-feet square with cropped corners. The 
control room and test section walls are provided with 
large windows for close viewing of the model. The 
tunnel is capable of operation at stagnation pressures 
from near vacuum to slightly above atmospheric and at 
Mach numbers from near zero up to about 1.2. Either 
air or a heavy gas (R-134a) can be used as the test me- 
dium. Both the density and test-section Mach number 
are continuously controllable. The present investigation 
was conducted in air near atmospheric conditions and at 
free-stream Mach numbers less that 0.30. 

Model 

The model used in the investigation is a modified 
version of a Froude-scaled semispan aeroelastic model 
of the V-22 tiltrotor that was used by Bell/Boeing to 
support the preliminary and full-scale design phases of 
the aircraft (Ref. 30). Upon completion of that series of 
tests, the Navy transferred the model to the Aeroelastic- 
ity Branch (AB) at NASA Langley under a loan agree- 
ment to be used as the experimental testbed of a tiltrotor 
aeroelastic research program that was initiated by AB in 
1994. The tiltrotor testbed has been designated the 
Wing and Rotor Aeroelastic Testing System (WRATS). 

General Characteristics 

The WRATS tiltrotor testbed as installed in the 
Langley Transonic Dynamics Tunnel tor this study is 
shown in figure 2. The model has a length scale factor 
of 1/5 and was designed to maintain full-scale values of 
the Froude, Lock, and Strouhal numbers when operated 
in air (Ref. 30). The wing and rotor are dynamically 
and aeroelastically scaled; the pylon is only dynamically 
scaled. The fuselage is rigid and only maintains the 
scaled external aerodynamic shape. The model is at- 
tached to a support structure that is effectively rigid, its 
lowest frequency being well above any important clastic 
mode frequency of the model. Simulation of the dis- 
tributed wing beamwise, chordwise, and torsional stiff- 


ness is provided by means of a hollow, rectangular cross 
section, composite spar having chordwise flanges. The 
4.6-foot spar, which lies along the calculated elastic axis 
of the wing, has segmented, nonstructural, aluminum 
aerodynamic fairings that provide the spanwise distribu- 
tion of airfoil contour. To provide surface continuity 
over the lifting surface of the wing, the space between 
the segments is filled with strips of foam rubber. The 
wing-tip-mounted pylon contains the transmission and 
gearbox components for the rotor drive system, the 
lower part of the mast, and the swashplate control sys- 
tem. Because these internal components can be treated 
as rigid, the pylon is scaled dynamically so as to pre- 
serve its overall mass and inertia properties. The pylon 
is attached to the wing tip by means of a “racetrack” 
spring assembly that simulates the combined stiffness of 
the full-scale conversion actuator and downstop lock 
mechanism. 

The 3-bladed, 7.6-foot diameter stiff-inplane ro- 
tor has a gimbalcd hub that is connected to the mast by a 
constant speed joint (2 coincident universal joints). The 
rotor yoke has 2.5 degrees of built-in precone and is 
flexible to allow the blade coning angle to adjust under 
centrifugal loading. The blades have a nonlinear distri- 
bution of built-in twist with an overall root-to-tip twist 
of 47.5 degrees (leading edge down). 

A large plywood panel through the vertical plane 
of symmetry of the model and attached to a support 
structure consisting of a spider-like arrangement of steel 
tubes was employed to seal the open back side of the 
semi-fuselage and to serve as a reflection plane. The 
structure also provided the offset needed to position the 
rotor axis near the centerline of the wind tunnel. 

The fundamental natural frequencies of the 
model (as measured on the test stand) with the pylon 
locked to the downstop in an airplane mode configura- 
tion are summarized in Table I. Typical wind-on natural 


Table I - System Natural Frequencies* 


Mode 

Frequency, Hz 

Wing Beam Bending 

5.83 

Wing Chord Bending 

8.67 

Wing Torsion 

12.0 ^ 

Pylon Yaw (2nd wing chord) 

19.4 

Blade 1st Elastic Flap 

7.20 

Blade 1st Elastic Lag 

12.5 

Blade Rigid Body Torsion 

95.6 

Blade 1st Elastic Torsion 

110. 


* Rotor on, flapping locked out, zero rpm. 
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frequencies of the rotor system (e.g., at V-lOO knots and 
Q. = 742 rpm) are: gimbal flapping at 0.85P, first cyclic 
lag at 1.2P, and collective flapping (coning mode) at 
about 1.7 P. The present test was run with the drive 
system disconnected so the collective lag mode had zero 
frequency. 

Active Swashplate Hardware 

Three high-frequency electro-hydraulic actuators 
having a bandwidth of 50 Hertz were used in the control 
system to drive the swashplate (fig. 4). The actuators 
are positioned azimuthally at the 12, 4, and 8 o’clock 
positions around and slightly aft of the nonrotating 
swashplate. These actuators were used to input the 
steady pilot commands as well as the swashplate inputs 
called for by the GPC active control algorithm. Three 
3000-psi servo-valves were used to meter hydraulic 
fluid to the actuators. These servo-valves were located 
inside the pylon near their corresponding actuators (see 
fig. 4). 



Figure 4. - Major hardware components of the GPC 
active control system. 


Feedback Signals and Instrumentation 

Feedback signals for the active control system 
were developed from the measured responses of the 
pylon/wing system. The suite of candidate responses 
included wing strain gages (beam, chord, and torsion) 
and six pylon accelerations (normal, lateral, and axial). 
The beam and chord strain gages are located 17.0 inches 
out from the root and the torsion strain gages are located 
23.0 inches out from the root. The wing bending and 
torsion gages were the sensors used for feedback. A 
number of other measurements were monitored and 
recorded either to trim the rotor or for flight safety, in- 


cluding gimbal flapping angles and downstop spring 
load. Pitch-link and blade loads were not measured. 

Control System Architecture 

A schematic diagram of the control system for 
active stability augmentation is shown in figure 5. The 
flapping response of the rotor is displayed 



Figure 5. - Schematic of active stability control system. 


to the pilot as IP lateral and longitudinal flapping of the 
rotor tip-path-plane. The pilot trims the windmilling 
rotor by setting both the collective pitch (to maintain the 
desired rotor speed) and the cyclic pitch (to maintain 
zero flapping with respect to the shaft). A mixer in the 
pilot control console calculates the commands to be sent 
to the individual actuators to trim the rotor based on the 
pilot stick commands. The GPC control system calcu- 
lates the actuator commands needed for stability aug- 
mentation. The steady OP (DC) swashplate actuator 
commands required for trim are summed with the oscil- 
latory actuator commands calculated by the GPC system 
(after those commands are converted from digital to 
analog signals). 

Data Acquisition and Processing 

The data acquisition system (DAS) that is part of 
the WRATS testbed is a 64-channel, PC-based system 
consisting of a National Instruments data acquisition 
board housed in a 750 MHz Pentium processor com- 
puter with a per-channel sampling rate of 1000 Hz. The 
WRATS DAS is controlled using a LabVIEW front 
panel that has been customized to support WRATS 
model testing. LabVIEW applications have been devel- 
oped for on-line FFT analyses, harmonic analyses, ex- 
citing the model via stick-stir inputs to the swashplate, 
damping estimations using moving block and logarith- 
mic decrement analyses, and plotting. The digitized 
data records are written in a binary file format to the 
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host PC’s hard disk for access during the test and to 
CDs for permanent storage. A MATLAB-based pro- 
gram is used for post-test data reduction and plotting. 


The GPC computer consisted of a 500 MHZ dual 
processor PC that included a dSPACE DSP1 103 add-in 
card for performing real-time data acquisition and com- 
putations. A sampling rate of 100 Hz was used for the 
GPC data channels. The GPC computational tasks were 
distributed among these three computing elements in the 
manner indicated in figure 6. The entire process was 
managed using ControlDesk, dSPACE’s graphical user 
interface (GUI), which was installed on CPU #2. The 
portion of the GPC algorithm that computes the control 
gain matrices d and p is written in MATLAB and runs 
on CPU #1 . The calculation of the control inputs u L that 
are sent to the swashplate actuators is made by executa- 
ble code that is installed on the DSP card. This code is 
generated by using Real-Time Workshop to convert a 
Simulink block diagram of the control law calculation to 
C code. The C code then goes to the dSPACE compiler, 
which generates the executable code that is placed on 
the DSP. The user specifies appropriate values for the 
parameters /, /;, h f „ h n w Vi and w r needed by the GPC 
algorithm and initializes the control gain matrices a 1 and 
P to zero using the ControlDesk front panel GUI that 
runs on CPU #2. The DSP card continuously collects 
the r-input and m-output data sets used by the GPC al- 
gorithm. On user command, the DSP sends the set of 
input/output data needed for system identification to 
CPU #1 where the SID computations are performed and 
the control gain matrices a and P are computed. The 
control gain matrices are automatically sent to the DSP, 
which uses the p latest data sets to (continuously) com- 
pute the control commands to be sent to the swashplate 
actuators. Re-identification of the system, if needed, is 
done closed-loop and on user command. 


CPU #2 


DSP Card < 


CPU#1 



Figure 6. - Distribution of computing tasks in GPC 
computer for active control studies. 


Test Conditions and System Configurations 

The GPC investigation constituted only a portion 
of a much broader aeroelastic investigation that was 
conducted with the WRATS model during the April 
2000 test in the TDT. One of the other objectives of 
that test was to establish new baseline stability bounda- 
ries for the model that included the hydraulic compo- 
nents that had been installed earlier for active controls 
testing. This baseline stability investigation was con- 
ducted with the pylon in both the on- and off-downstop 
configurations, over a range of rotor rpm, in both air and 
heavy gas. 

All GPC testing was conducted with the pylon 
oriented in the airplane mode of flight and locked to the 
wing downstop. The rotor was unpowered (wind- 
milling). The stability of the model with the GPC-based 
algorithm was measured over a range of tunnel air- 
speeds up to a maximum of 120 kt, for a single rotor 
rotational speed (888 rpm), and two values of kinematic 
pitch-flap coupling (5 3 ) angles (-15° and - 45°). Be- 
cause this was the first wind-tunnel test of the GPC al- 
gorithm, a conservative approach was adopted wherein 
all closed-loop testing was conducted at conditions that 
were slightly within (more stable than) the correspond- 
ing open-loop stability boundaries. 

Test Procedure 

The primary pylon/wing structural frequencies of 
interest lie between 4 and 25 Hz, with the critical wing 
beam mode frequency varying from 5.8 Hz down to as 
low as 4.5 Hz with increasing airspeed. Based on the 
implementation guidelines discussed earlier, the Nyquist 
frequency f N was taken as 30 Hz, the sampling fre- 
quency f s was set to 1 00 Hz, and the number of time 
steps / of input/output data used for SID was taken as 
300. 


For this test, the control horizon h c was set equal 
to the prediction horizon h r , h j} was set equal to /?, the 
order of the ARX model, and the response weight w r 
was set to unity. The assumed order of the ARX model 
(p) and the control weight (w\.) were varied to tune the 
controller to the nuances of the WRATS model in the 
on-downstop configuration with 63 = — 15° while operat- 
ing at an airspeed of 80 kt, which is well below the 
open-loop flutter speed of that configuration. The tun- 
ing process was conducted at both 770 and 888 rpm and 
led to the selection of the following best values for these 
quantities: p = 30 and w c = 1 .0. Active control was in- 
troduced into the fixed-system swashplate using three 
hydraulic servo-actuators mounted aft of the swashplate 
inside the pylon fairing. Wing bending (vertical and 
ehordwise) and torsion strain gages located near the 
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wing root were the sensors used to provide the response 
measurements needed by the active control algorithm. 
Thus, there were three control inputs (r = 3) and three 
feedback outputs (m = 3). 

In addition to its use to actively control the rotor 
for GPC stability augmentation, the active swashplate 
was also used to impose “stick stirs 5 ' to excite the wing 
modes of interest for damping determination. These 
stick-stir commands were superimposed on the GPC 
and pilot control commands sent to the swashplate ac- 
tuators. Control of the stick stir is integrated into the 
GUI that is part of the WRATS DAS software. The 
procedure employed for the GPC testing is as follows: 
(1) With GPC off, the tunnel airspeed is brought up to a 
selected value; (2) GPC is turned on (the associated 
reduction in model response is visibly quite dramatic); 
(3) Data acquisition is initiated; (4) A stick-stir excita- 
tion is imposed at the frequency of the wing mode of 
interest; (5) The resulting wing beam (or other) response 
is observed visually by the test engineer and the stir 
terminated when sufficient amplitude has been attained; 
( 6 ) The response is allowed to decay freely to its steady- 
state value ; and (7) Data acquisition is turned off. The 
total data acquisition time is typically 5-8 seconds long. 
The free-decay portion of the recorded time history is 
then analyzed by moving block and log decrement 
methods to obtain an estimate for the modal damping 
and the damped natural frequency. This process is re- 
peated until five acceptable estimates of the damping 
are obtained. GPC is then turned off and the entire pro- 
cedure beginning with step ( 1 ) is repeated for all of the 
airspeeds of interest. 

The performance and robustness of the GPC al- 
gorithm were assessed in transient conditions for S 3 = 
-15° by rapidly changing tunnel airspeed and rotor rota- 
tional speed, both singly and in combination. The rotor 
speed was varied between 770 and 888 rpm and the 
tunnel airspeed between 60 and 145 kt. 

RESULTS AND DISCUSSION 

The GPC-based active control system was found 
highly effective in increasing the stability (damping) of 
the critical wing mode for all of the configurations of 
the model tested. In particular, GPC was able to consis- 
tently yield a closed-loop system in which the critical 
wing mode had a minimum of about four-percent modal 
damping over the range of test conditions investigated, 
without visible degradation of the damping in the other 
modes. The algorithm was robust with respect to its 
performance in the tracking of rapid changes in both the 
rotor speed and the tunnel airspeed. System identifica- 
tion done at a low-speed flight condition was generally 


sufficient for a wide range of rotor speeds and tunnel 
velocities. 

An indication of the effectiveness of the GPC- 
based active control system in increasing the stability 
(damping) of the critical wing beam mode of the model 
is given in figures 7 and 8 . The figures show a com- 
parison of the measured open-loop and closed-loop 
wing beam mode damping versus airspeed for 83 values 
of -45° and -15°, respectively. The mean and standard 



Figure 7. - Comparison of wing beam mode damping 
versus airspeed with GPC on and off for 5^ = —45°. 



Figure 8 . - Comparison of wing beam mode damping 
versus airspeed with GPC on and off for 8^ = — 15°. 

deviation of the five values of damping that were meas- 
ured at each airspeed were determined. The mean val- 
ues of damping are indicated by the symbols, and the 
vertical lines indicate the standard deviations about 
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these mean values. The plotted curves are the result of a 
spline fit applied to the mean values. As can be seen, 
with GPC turned on the wing beam mode damping is 
considerably higher than with GPC off over the entire 
range of tunnel airspeeds tested. It should be remarked 
that with GPC on, it was often difficult to estimate the 
damping from the free-decay responses after termina- 
tion of the stick-stir because the large damping levels 
resulted in very few cycles of motion (sometimes as few 
as one or two) being available for the damping calcula- 
tions. This is the reason for the large standard deviation 
in the measured values of the closed-loop damping. In 
contrast, because the open-loop damping of the wing 
beam mode is small, the standard deviation of the meas- 
ured damping is correspondingly very small. It should 
also be remarked that the effectiveness of GPC was 
such that it tried to quell the imposed stick-stir excita- 
tion, thus requiring considerable amplitude of excita- 
tion. The reason for the drop-off and then rise in 
closed-loop damping at 100 knots in figure 8 is not 
known. 

The measured time histories of the wing beam 
bending moments during and after a typical stick-stir 
excitation with GPC on and off are shown in figures 9 
and 10. Figure 9 corresponds to the 75-knot airspeed of 
figure 7 and figure 10 corresponds to the 120 -knot air- 
speed of figure 8 . The closed-loop time histories clearly 
illustrate the effectiveness of the GPC-based active con- 
trol system to rapidly reduce the response. 



Figure 9. - Effect of GPC on damped response of wing 
beam bending mode at 75 kt for 83 = -45°. 



1 1 ’ \ l 1 1 1 1 1 1 

0 1.5 3.0 4.5 


Time, sec 

Figure 10. - Effect of GPC on damped response of wing 
beam bending mode at 120 kt for 83 = -15°. 

Another example of the effectiveness of GPC is 
given in figure 1 1 , which shows a time history of the 
wing beam bending moment while repeatedly cycling 
GPC on and off for the 120-knot airspeed case of figure 
8 . The response is seen to adapt quickly to the changing 
conditions, suggesting a certain level of robustness. 



Figure II.- Time history of wing beam bending mo- 
ment while cycling GPC on and off at 120 kt for 83 = 
-15°. 

An indication of the swashplate pitch angles that 
are associated with the GPC commanded actuator inputs 
is given in figure 12. The figure shows the time history 
of the longitudinal and lateral cyclic pitch angles for the 
120 -knot airspeed condition of figure 8 from a time near 
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the end of a stick-stir excitation to a time well into the 3) All system identification and control law 

steady-state condition. The steady-state oscillatory cy- computations were done on-line, permitting rapid adap- 
clic pitch angles are seen to be modest. tation to changing conditions. 



Time, sec 


4) The swashplate angles required were modest, 
generally less than 0.3 degrees. 

These results, in combination with equally suc- 
cessful results obtained in an earlier ground resonance 
test, suggest that a GPC-based active control system is a 
viable candidate for stability augmentation in advanced 
tiltrotor systems. Additional wind-tunnel tests of the 
WRATS model are planned to evaluate the GPC meth- 
odology over a broader range of operating conditions 
for the baseline stiff-inplane gimbaled rotor, as well as 
an advanced soft-inplane design. 
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